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K^ , Abstract 

^ Cy , More than half a milhon minimum-bias events of LHC coUision data 

were collected by the ATLAS experiment in December 2009 at centre-of- 
mass energies of 0.9 TeV and 2.36 TeV. This paper reports on studies of the 
OO . initial performance of the ATLAS detector from these data. Comparisons 

between data and Monte Carlo predictions are shown for distributions of 
several track- and calorimeter-based quantities. The good performance 
^ ■ of the ATLAS detector in these first data gives confidence for successful 

O ' running at higher energies. 
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1 Introduction 



In December 2009, the ATLAS detector P recorded data from a first series 
of LHC runs at centre-of-mass energies of 0.9 TeV and 2.36 TeV. When the 
'.^ ■ beams were colliding and declared to be stable by the LHC operators, all main 

l/^ I detector components were fully operational and all levels of the trigger and 

CN ' data acquisition system performed as expected, assuring smooth and well mon- 

"^ I itored data taking. The data sample at 0.9 TeV contains nearly 400 000 events 

l/^ ' recorded with high-quality calorimeter and tracking information, corresponding 

^D I to an integrated luminosity of approximately 9 |xb~^ [2]. The data at 2.36 TeV, 

36 000 events, which are only used here for calorimeter studies, correspond to 
approximately 0.7 |J.b~^. These data sets do not contain very many high-px 
objects, and therefore do not correspond to the environment for which ATLAS 
was designed. 
5J] ' The ATLAS detector was thoroughly commissioned and initial calibration 

5^ I and performance studies were done using cosmic ray data recorded during 2008 

and 2009. Performance close to design goals was obtained for the different 
detector components, details can be found in Refs. [3J|11[S1|B]. 

This paper presents performance established with data taken in first colli- 
sions in 2009. The detector components are outlined in Section [21 The trigger 
and data acquisition performance together with the initial event selection are 
discussed in Section [31 the simulation to which the data are compared is ex- 
plained in Section [31 The performance of the inner tracking system is reviewed 
in Section [5l and the combined analysis of calorimeter data and tracking infor- 
mation to study electrons and photons is discussed in Section [SI Studies of jets 
and missing transverse energy, E!f^^^, using the calorimeter cells are presented 



in Sections [7] and |H1 Finally, kinematic distributions of the first reconstructed 
muon candidates are shown in Section [3] 

2 The ATLAS Detector 

The ATLAS detector [Tl covers almost the entire solid angle around the nominal 
interaction point and comprises the following sub-components: 

• An inner tracking system: operating inside an axial magnetic field of 
2 T, it is based on three types of tracking devices. These are an outer 
tracker using straw tubes with particle identification capabilities based on 
transition radiation (Transition Radiation Tracker, TRT), a silicon strip 
detector (Semiconductor Tracker, SCT) and an innermost silicon pixel 
detector (Pixel). 

• A hybrid calorimeter system: for the electromagnetic portion (EM), the 
hadronic end-cap (HEC) and the forward calorimeter (FCal) a liquid ar- 
gon (LAr) technology with different types of absorber materials is used. 
The central hadronic calorimeter (Tile) is a sampling calorimeter with 
steel as the absorber material and scintillator as the active medium. The 
electromagnetic sections use an accordion geometry to ensure fast and 
uniform response. A presampler detector, to correct for energy losses in 
the upstream material, is installed in front of the EM calorimeter in the 
range \r]\ < 1.80 

• A large muon spectrometer: an air-core toroid system generates an average 
field of 0.5 T (1 T), in the barrel (end-cap) region of this spectrometer, 
resulting in a bending power between 2.0 and 7.5 Tm. Over most of 
the 77-range, tracks are measured by Monitored Drift Tubes (MDT); in 
the high ?7-regime the closest of four wheels to the interaction region is 
instrumented with Cathode Strip Chambers (CSC). Trigger information 
is provided by Thin Gap Chambers (TGC) in the end-cap and Resistive 
Plate Chambers (RPC) in the barrel. 

• Specialized detectors in the forward region: two dedicated forward de- 
tectors, the LUCID Cherenkov counter and the Zero Degree Calorime- 
ter (ZDC). In addition the BPTX, an electrostatic beam-pickup which 
monitors the timing of the beam near ATLAS and two scintillator wheels 
(MBTS) were mounted in front of the electromagnetic end-caps to provide 
trigger signals with minimum bias. 

3 Data- Taking Performance and Event Selection 

The ATLAS operating procedure in 2009 maintained the calorimeters and TRT 
in standard operating conditions, but the silicon trackers and muon chambers 



^The origin of the coordinate system used to describe the ATLAS detector is the nominal 
interaction point. The positive x axis is defined as pointing from the interaction point to 
the centre of the LHC ring, the positive y axis is defined as pointing upwards and the beam 
direction defines the z axis of a right-handed coordinate system. Transverse momenta are 
measured in the x-y plane with radius r. Polar {9) and azimuthal (rp) angles are measured 
with respect to this reference system. The pseudorapidity is defined as r; = — lntan(6/2). 



Table 1: Luminosity- weighted fraction of the time during stable beam operation 
for which the different detectors were able to take data under nominal conditions. 

Pixel SCT TRT LAr Tile MPT RFC TGC 
Efficiency [%] 80.9 86.2 100 99.0 100 87.4 88.6 84.4 



were at a reduced or 'standby' voltage until after stable beams were declared by 
the LHC. Most of the studies in this paper required the tracking detectors to be 
in operating conditions and used approximately 400 000 events while the E^^^^ 
studies only required the calorimeters and used some 600 000 and 36 000 events 
at 0.9 TeV and 2.36 TeV respectively. The luminosity-weighted availability of 
the various sub-detectors (see Section 13. 3p during stable beam operations is 
summarized in Table [1] 

3.1 Trigger/DAQ System 

The ATLAS trigger and data acquisition (TDAQ) is a multi-level system with 
buffering at all levels [TI. Trigger decisions are based on calculations done at 
three consecutive trigger levels. While decisions at the first two levels are pend- 
ing, the data acquisition system buffers the event data from each sub-detector. 
Complete events are built after the second level decision. The first level trigger 
(LI) is largely based on custom built electronics. It incorporates timing from 
the BPTX and coarse detector information from the muon trigger chambers 
and the trigger towers of the calorimeters, along with multiplicity information 
from the MBTS scintillators and the ATLAS forward detectors, LUCID and 
ZDC. The LI system is designed to select events at a rate not exceeding 75 kHz 
from an input rate of 40 MHz and identify regions-of- interest (Rols), needed by 
the high-level trigger system (HLT), for potentially interesting physics objects. 
HLT runs on a processor farm and comprises the second level (L2) and the third 
level trigger or Event Filter (EF). The L2 system evaluates event characteris- 
tics by examining the Rols using more detector information and more complete 
algorithms. The EF analyzes the L2 selected events again looking at the Rols' 
measurements. The analysis of the complete event is also possible. The output 
rate is reduced to approximately 200 Hz. 

All trigger decisions used were made by the LI systems as the luminosity 
was low. Minimum-bias events were triggered by a coincidence between the the 
BPTX and signals indicating hits in one or both of the MBTS scintillator wheels. 
However, the functionality of the main L2 and EF algorithms, for example those 
used for inner track reconstruction and jet finding, were validated by running 
in 'passthrough' mode, i.e. calculating relevant L2 and EF decisions without 
rejecting events. Important distributions like the vertex position (Fig. [TJa)) 
were monitored online and different data streams for fast analysis feedback based 
on LI trigger decisions were provided. For the December 2009 running period 
the data-taking efficiency of the overall Trigger/DAQ system, as illustrated in 
Fig. mb), averaged to 90%. 
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Figure 1: (a) Reconstructed z-vertex distribution calculated online by the higher 
level trigger for monitoring purposes. The width includes a small contribution 
from the experimental resolution, (b) Data-taking efficiency for periods with 
two circulating beams in December 2009. 
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Figure 2: Time difference AtuBTS of hits recorded by the two MBTS scintillator 
wheels mounted in front of the electromagnetic end-cap wheels on both sides of 
the ATLAS detector; (a) time difference without any selection and (b) requiring 
a well-reconstructed vertex. 



3.2 Event Selection 

To select collision candidates and remove beam-related background two different 
strategies were employed: 

• For those studies based mainly on track information, the presence of a 
primary vertex, reconstructed using at least three tracks with sufficient 
transverse momenta, typically pr > 150 MeV, and a transverse distance 
of closest approach compatible with the nominal interaction point are 
required. This selection strategy, which was used in Ref. [2 , uses events 
triggered by a single hit in one of the two MBTS scintillator wheels. 

• Alternatively, the selection is based on the timing difference of signals de- 
tected on both sides of the ATLAS detector. Coincident signals, within 
a time window of 5 or 10 ns from either the electromagnetic calorime- 
ters (end-cap or FCal) or from the two MBTS wheels, respectively, are 
required. The event must again be triggered by an MBTS signal. In case 
no timing coincidence is found, a two hit MBTS trigger with at least one 
hit per side is required. 

The detailed track quality criteria used for the first strategy vary slightly for 
the different studies presented and are described later when appropriate. 

Without any event selection the MBTS-triggered events contain backgrounds 
from beam-related events as shown in Fig. [Ifa), where the time difference, 
AiMBTS 7 of MBTS signals recorded on both sides of the ATLAS detector is de- 
picted. For events coming from the interaction point AtMBTS is small. Beam- 
related background produced upstream or downstream should have A<mbts 
around 25 ns, with the sign giving the direction. Eighty percent of single beam 
events are missing timing information on one or both sides and are therefore not 
shown. Requiring a well-reconstructed vertex with track quality requirements 
reduces the beam-related background by more than three orders of magnitude 
while retaining genuine collision events (Fig. HJb)). There are twelve single- 
beam events which meet this vertex requirement, but all of them are missing 
timing information and are not shown. 

3.3 Luminosity Measurement 

The luminosity during the 2009 ATLAS data-taking period was estimated offline 
based upon the timing distributions measured by the MBTS. Events with signals 
detected on opposite ends of the ATLAS detector in the MBTS are counted. 
After background subtraction the luminosity is calculated using the number 
of events with a timing difference consistent with particles originating from 
the interaction point (see Fig. [5]), the expected minimum-bias cross section 
and the event selection efficiency determined from data and Monte Carlo. The 
MBTS detector is used for the absolute luminosity determination because of 
its high trigger efficiency for non-diffractive events 0. The uncertainty on the 
luminosity, dominated by the understanding of the modelling of inelastic pp 
interactions, is estimated to be around 20%. 

Figure [3] shows the luminosity as a function of time, as measured using the 
MBTS system as well as with three other techniques: timing in the LAr, the 
LUCID relative luminosity monitor and particle vertices reconstructed online. 
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Figure 3: Instantaneous luminosity measured by the MBTS and LAr, with 
superimposed the LUCID and HLT vertex counting estimates normahzed in 
such a way to give the same integrated luminosity as measured with the MBTS 
system. All measurements are corrected for TDAQ dead-time, except LUCID 
which is free from dead time effects. The short luminosity drop at 14:15 is due 
to inhibiting the trigger for ramping up the silicon detectors after declaration 
of stable LHC beams. 



The LAr technique has a slightly larger systematic uncertainty in the accepted 
cross section than that from the MBTS and produces a result which agrees to 
4%. The other methods are normalized to the MBTS measurement. 



4 Monte Carlo Simulation 

Monte Carlo samples produced with the PYTHIA 6.4.21 j7] event generator are 
used for comparison with the data. ATLAS selected an optimized parameter 
set [S], using the px-ordered parton shower, tuned to describe the underlying 
event and minimum bias data from Tevatron measurements at 0.63 TeV and 
1.8 TeV. The parton content of the proton is parameterized by the MRST LO* 
parton distribution functions 0. 

Various samples of Monte Carlo events were generated for single-diffractive, 
double-diffractive and non-diffractive processes in pp collisions. The differ- 
ent contributions in the generated samples were mixed according to the cross- 
sections calculated by the generator. There was no contribution from cosmic 
ray events in this simulation. All the events were processed through the ATLAS 
detector simulation program [10 , which is based on GEANT4 11 . This sim- 
ulation software has also been systematically compared to test-beam data over 
the past decade (see e.g. Ref. [12) and it was constantly improved to describe 
these data. After the detector simulation the events were reconstructed and 
analyzed by the same software chain also used for data. 

The beam position and size, which did not correspond precisely to those used 
in the simulation prepared beforehand, have a significant impact on some distri- 



butions, particularly for detectors close to the interaction region. The length of 
the luminous region, as seen in Fig. [IJa), is approximately half that expected, 
and the simulated events were re-weighted to match this. The transverse offset 
in the simulation of about 2 mm cannot be corrected for by this method. 

The distributions presented in this paper always show the simulated sample 
normalized to the number of data events in the figure. 

5 Tracking Performance 

The inner tracking system measures charged particle tracks at all (jj and with 
pseudorapidity \ri\ < 2.5. The pixel detector is closest to the beam, covering 
radial distances of 50-150 mm with three layers both in the barrel region and 
in each end-cap. The innermost Pixel layer (known as the B-layer) is located 
just outside the beam pipe at a radius of 50 mm. The pixels are followed, 
at radii between 299 - 560 mm, by the silicon strip detector known as the 
SCT. This provides 4 (barrel) or 9 (end-cap) double layers of detectors. The 
Pixels are followed, for radii between 563 - 1066 mm, by the TRT. The TRT 
straw layout is designed so that charged particles with transverse momentum 
Pt > 0.5 GeV and with pseudorapidity \r]\ < 2.0 cross typically more than 30 
straws. The intrinsic position resolutions in r<j> for the Pixels, the SCT and the 
TRT are 10, 17 and 130 |J.m, respectively. For the Pixels and the SCT the other 
space coordinate is measured with 115 and 580 [im accuracy, where the SCT 
measurement derives from a 40 mrad stereo angle between the two wafers in a 
layer. 

5.1 Hits on Tracks 

The sample of minimum-bias events provides approximately two million charged 
particles with px over 500 MeV through the central detectors of ATLAS. Their 
trajectories in the inner detector were reconstructed using a pattern recognition 
algorithm that starts with the silicon information and adds TRT hits. This 
'inside-out' tracking procedure selects track candidates with transverse momenta 
above 500 MeV [IS]. Two further pattern recognition steps were run, each 
looking only at hits not previously used: one starts from the TRT and works 
inwards adding silicon hits as it progresses and the other repeats the first step, 
but with parameters adjusted to allow particle transverse momenta down to 
100 MeV. The multiple algorithms are necessary partly because a 100 MeV 
Pt charged particle has a radius of curvature of about 17 cm in the ATLAS 
magnetic field and will not reach the TRT. 

The track selection requirements vary slightly among the analyses presented 
here. A typical set of selections is that charged particle tracks are required 
to have pt>0.5 GeV, >1 Pixel hit, >6 SCT hits and impact parameters with 
respect to the primary vertex of | do I < 1-5 nun a-nd |zosin0| < 1.5 mm. The 
transverse impact parameter, do, of a track is its distance from the primary 
vertex at the point of closest approach when projecting into the transverse 
plane, signed negative if the extrapolation inwards has the primary vertex to 
the right, zq is the longitudinal distance at that point. 

The hit distributions in the silicon detectors as a function of (f) are shown 
in Fig. m for tracks passing these requirements in data and simulation. The 
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Figure 4: A comparison of data and simulation in the average number of hits 
in (a) the Pixels and (b) the SCT versus (f> on selected tracks. Comparable 
distributions versus rj can be found in Ref. [2]. 

fluctuations seen in correspond to non-responsive detector modules which are 
modelled in the simulation. A small mis-match between data and simulation 
arises because the simulated beam had a transverse displacement of about 2 mm 
from the true position, as discussed in Sectional 

The efficiency of the individual TRT straws is displayed in Fig.[5]as a function 
of the distance of the test track to the wire in the centre of the straw. The 
efficiency for data and simulation, barrel and end-cap, has a plateau close to 
94%. 
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Figure 5: TRT hit efficiency as a function of the distance of the track from the 
wire in the centre of the straw in (a) the barrel and (b) the end-caps. 



The alignment of the tracking detectors benefited from the precision con- 
struction and survey followed by an extended period of data taking using cosmic 
ray muons [5]. The alignment was improved using the 0.9 TeV collision data, 
although the particles have rather low momentum and therefore their tracks 
suffer from multiple scattering. The quality of the alignment can be checked by 
the study of the residuals, which are defined as the measured hit position minus 
that expected from the track extrapolation. 

Unbiased residuals between tracks and barrel TRT hits are plotted in Fig. [HI 
This figure is made using charged particles with pT > f GeV with over 6 hits in 
the SCT and at least f 4 in the TRT. The equivalent Gaussian width is extracted 
from the full- width at half maximum. The end-cap shows a resolution somewhat 
worse than simulation, while in the barrel part of the detector, where a higher 
cosmic ray flux can be used for alignment, data and simulation are in close 



agreement. 
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Figure 6: Unbiased residual distributions in the TRT barrel (a) and end-caps 
(b). The data points are in filled circles and the simulation in empty ones. 

Unbiased x residuals from the silicon detectors are shown in Fig. [71 where 
X refers to the more precise local coordinate on the detector. Charged particles 
are selected to have pt> 2 GeV. The equivalent Gaussian width is extracted 
from the full-width at half maximum. The width of the resulting distributions 
in data are within about 15% of those found in a simulation with no alignment 
errors, showing that the remaining impact on the residual widths from imperfect 
alignment in data is at the level of approximately 10-15 |J.m for the pixels and 
of 20 \iia for the SCT. 
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Figure 7: The distributions of the silicon detector unbiased residuals for (a) the 
pixel barrel, (b) the pixel end-cap, (c) the SCT barrel, (d) the SCT end-cap. 
The data arc in solid circles, the simulation, which has a perfect alignment, is 
shown with open ones. 
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Figure 8: The K^ candidate mass distribution using impact parameter and life- 
time selections. The simulated signal and background are separately normalized 
to the data. 

The momentum scale and resolution of the tracker, and energy loss with in, 
were all investigated by studying the Kg to tt+tt" decay. The reconstruction 
requires pairs of oppositely-charged particles compatible with coming from a 
common vertex. This vertex, in the transverse plane, must be more than 0.2 mm 
from the primary vertex. The cosine of the angle between the flight path relative 
to the primary vertex and the momentum vector of the candidate, cos6k, is 
required to exceed 0.8. The invariant mass distribution, calculated assuming 
that both charged particles are pious is shown in Fig. |8l The simulated signal 
and background are separately normalized to the data, and the position and 
width of the Kg mass peak are fitted using a Gaussian. The peak in data is at 
771^^ — 497.5 ± 0.1 MeV, in agreement with the PDG average ^^. 

In order to test the momentum scale and resolution of the detector the 
reconstructed pions in the simulation are adjusted by parameters ^tr, which 
scales the 1/pt, and atr, a Gaussian smearing on /utr- The values of these 
parameters which best fit the observed Kg mass and width in the barrel region 
are /Xtr = 1.0004 ±0.0002 and a^ = 0.0040 ±0.0015. Thus the momentum scale 
for these barrel charged particles is known at better than the one per mille level, 
which is as expected from the accuracy of the solenoid magnet field-mapping 
performed before installation of the inner detector [15]. This, and subsequent 
Kg studies, use a tighter cut of 0.99 on cos 9k- 

In the end-cap regions there is evidence for a degraded resolution, especially 
at low momentum. Charged particles with px below 500 MeV require a crtr of 
0.024±0.004 and 0.022±0.004 in the negative and positive end-caps, respectively, 
to match the data, suggesting some material is missing in the description of the 
end-caps. The momentum scale in the end-caps is compatible with the nominal 
within errors of 1 to 2 per mille. 
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The Kg peak was also used to investigate the amount of material in the 
inner tracker as a function of radius. The mass reconstructed in data, divided 
by that found in simulation, is shown in Fig. IHlas a function of decay radius. 
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Figure 9: The fitted A'g mass divided by the value found in nominal MC simu- 
lation as a function of the reconstructed decay position. The filled circles show 
the data, and the open symbols are for simulation samples with approximately 
10% and 20% more silicon tracker material added. The horizontal dotted line 
is to guide the eye. 

Deviations of this ratio from unity would expose differences between the real 
detector and the model used for simulation. The results for special simulation 
samples with approximately 10% and 20% fractional increase in the radiation 
length of the silicon systems included by increasing the density of some of the 
support structures are also shown in Fig. [SI These results suggest that discrep- 
ancies of material between the data and the simulation must be significantly 
smaller than 10% of the material thickness in the inner silicon barrels. 

5.3 dE/dx and 0(1020) Identification 

One feature of the Pixel tracking system is a time-over-threshold measurement 
for the signal which was used to extract the specific energy loss dE/dx. Tracks 
with more than one Pixel hit were studied and the mean dE/dx was found for 
each after the highest was removed to reduce the effect of Landau fluctuations. 
Figure [TU] shows the distribution observed in the data. Bands corresponding to 
different particle species are clearly visible. 

In the observation of — >■ K^K~ identification of the K^ reduces the com- 
binatorial background. The identification of kaon candidates through dE/dx 
proceeds by finding the probability density functions of pions (ppion), kaons 
(Pkaon) and protons (pproton) in the simulation as a function of momentum and 
dE/dx. This is done via fitting the observed value in simulation using a Gaus- 
sian function whose parameters are momentum dependent. The simulation 
models the data with an accuracy of about 10%. 

The tracks used in the reconstruction of the </> meson must have more than 
one hit in the Pixel system and an impact parameter within Scr of the primary 
vertex. The track fit was re-run using the kaon mass hypothesis for the energy 
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Figure 10: The dE/dx measured in data as a function of momentum. 



loss. The simulation shows that after re-fitting the kaon momenta are underesti- 
mated by up to 10 MeV and a corresponding correction is applied. This changes 
the reconstructed </) mass by approximately 0.3 MeV. All oppositely charged par- 
ticle pairs where both momenta, reconstructed under the kaon hypothesis, are 
below 800 MeV are considered. 
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Figure 11: The measured and simulated mass spectra of K^K^ pairs. The 
4> peak is fitted with a Breit-Wigner with a fixed width convoluted with a 
Gaussian. Both kaons must be identified through the dE/dx measurement. 

Figure [TT] shows the resulting mass distribution for the K~^K~ candidate 
pairs, selected using charged particles with 200 < px < 800 MeV and a kaon 
dE/dx tag. The selection cuts were chosen to yield optimal signal significance 
on simulated events; a measure which was greatly improved using the dE/dx 
information. 
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The background and signal levels in the simulation were scaled independently 
to match the data. The fit allowed the mass and experimental resolution to vary, 
while keeping the natural width fixed to the PDG JJ] average. The mass was 
found to be 1019. 5±0. 3 MeV, in agreement with the expected value. The fitted 
experimental resolution in data is 2.5±0.5 MeV and matches the 2.4±0.3 MeV 
found in Monte Carlo simulation. 

5.4 Secondary Vertex Tagging 

An important role of the tracking system is the identification of heavy flavour 
hadrons. There are several tagging algorithms developed in ATLAS. Some per- 
formance figures for two algorithms, the impact parameter and the secondary 
vertex tagging algorithm, are presented in the following. 

The transverse impact parameter, do, is a key variable for discriminating 
tracks originating from displaced vertices from those originating from the pri- 
mary vertex. For studies of track impact parameters the do was calculated with 
respect to a primary vertex which was fitted excluding that track in order to 
remove bias. 

In order to study the effect of material on the do resolution. Fig. [T2T a) shows 
a^{do) versus 1 / {p-^ sin'^ 9) for data and simulation using all selected charged 
particle tracks. The quantity a{do) is determined by fitting the do distribution 
in each bin of l/(p^ sin'^ 9) with a Gaussian within ±2a{do) about its mean. The 
data lie approximately on a straight line, as is expected if the scattering material 
is on a cylinder and the match of the slope with the simulation implies a good 
description of the material of the inner detector. It should be noted that the 
intercept on the y axis has a contribution from the primary vertex resolution. 
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Figure 12: (a) The variance of the do distribution as a function of l/{p^ sin'^ 9) 
of the tracks for data (solid points) compared to the nominal simulation (open 
points). A straight line fit to the data points is also shown, (b) The lifetime- 
signed impact parameter significance. 

The track selection for the 5-tagging algorithms is designed to select well- 
measured particles and reject badly measured tracks, tracks from long-lived 
particles (Xg, A and other hyperon decays), and particles arising from material 
interactions such as photon conversions or hadronic interactions. 
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Table 2: Track selection criteria used for the impact parameter and secondary 
vertex tagging algorithms. 



Impact parameter Standard vertex Loose vertex 



Pt 


> 1 GeV 


do 


< 1 mm 


zq sin 9 


< 1.5 mm 



> 0.5 GeV 

< 2 mm 

< 2 mm 



> 0.5 GeV 

< 10 mm 

< 50 mm 



The track selection used by the impact parameter tagging algorithm is sum- 
marized in the first column of Table [5] Slightly different selections are used by 
the secondary vertex algorithm (second column of Tabled]). 

Calorimeter jets (see Section [T]) are the reconstructed objects the tagging 
algorithms are typically applied to. Their direction is taken as estimator of the 
putative heavy flavour hadron direction. The impact parameter is then signed 
by whether the track perigee, relative to the jet direction, suggests a positive or 
negative flight distance. The distribution of the lifetime-signed impact param- 
eter significance for tracks in jets is shown in Fig. I12f b). 

Reconstructing explicitly the secondary decay vertices of heavy flavour had- 
rons adds substantial tagging information. There are expected to be few b- 
quarks which can be tagged in this data set, so the algorithm was run with the 
standard as well as loose settings, as described in the second and third columns 
of Table [21 respectively. The loose setting selects vertices originating from K^ 
as well as from b-hadron decays whereas in the standard configuration any pair 
of tracks consistent with a K^, A or photon conversion is explicitly removed. 

Secondary vertices are reconstructed in an inclusive way starting from two- 
track vertices which are merged into a common vertex. Tracks giving large x^ 
contributions are then iteratively removed until the reconstructed vertex fulfils 
certain quality criteria. 

The mass distribution of the resulting vertices for the loose configuration, 
assuming a pion mass for each track, is shown in Fig. 1131 
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Figure 13: The vertex mass distribution for all secondary vertices with pos- 
itive decay length selected in data. The expectation from simulated events, 
normalized to the number of jets in the data, is superimposed. 
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Running the algorithm in the standard configuration results in the recon- 
struction of 9 secondary vertices with positive decay length significance. This 
is in good agreement with the 8.9 ± 0.5(stat.) vertices expected from the same 
number of jets, 10 503, in non-diffractive minimum-bias simulation. The vertices 
reconstructed with the standard version of the tagging algorithm are predom- 
inantly those with higher masses as the low-mass region is dominated by Kg 
mesons. 



ATLAS 

EXPERIMENT 




Figure 14: An event containing a secondary vertex selected by the secondary 
vertex algorithm. The pixel detector can be seen on the left and an expansion 
of the vertex region on the right. Unassociated hits, in a lighter colour, are 
predominantly due to unreconstructed particles such as those with transverse 
momenta below 0.5 GeV. 

An event display of the highest-mass candidate is shown in Fig. [141 The 
secondary vertex consists of five tracks and has a mass of 2.5 GeV. The vertex 
is significantly displaced from the primary vertex, with a signed decay length 
significance L/ a{L) — 22. From the vertex mass, momentum and L a proper 
lifetime of 3.1 ps is estimated. The data was also tested by the impact-parameter 
based b-tagging algorithm and this jet is assigned a probability below 10"'' for 
originating from a light quark jet. 

5.5 Particle Identification using Transition Radiation 

The TRT provides substantial discrimination between electrons and pions over 
the wide energy range between 1 and 200 GeV by utilizing transition radiation 
in foils and fibres. The readout discriminates at two thresholds, the lower set 
to register minimum-ionising particles and the higher intended for transition 
radiation (TR) photon interactions. The fraction of high-threshold TR hits as a 
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function of the relativistic 7 factor is shown in Fig. [12] for particles in the forward 
region. This region is displayed because there are more conversion candidates 
and they have higher momenta than in the barrel. 
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Figure 15: The fraction of high-threshold transition radiation hits on tracks as 

a function of the relativistic 7 factor (see text for details). 

The high-7 part of the distribution is constructed using electrons from pho- 
ton conversions while the low-7 component is made using charged particle tracks 
with a hit in the B-layer and treating them as pious. All tracks are required to 
have at least 20 hits in the TRT. The photon conversions are found similarly 
to those in Section 15771 with at least one silicon hit, but the transition radiation 
electron identification was not applied to the electron that was being plotted. 
To ensure high purity (about 98%), the conversion candidates are also required 
to have a vertex more than 40 mm away from the beam axis. The pion sample 
excludes any photon conversion candidate tracks. 

5.6 Tracking EfRciency for Level-2 Trigger 

The L2 track trigger is one component of the HLT whose performance can be 
tested with current data. The trigger runs custom track reconstruction algo- 
rithms at L2, designed to produce fast and efficient tracking using all tracking 
subdetectors. Tracking information forms an integral part of many ATLAS 
triggers including electron, muon and tan signatures |16| . These use LI infor- 
mation to specify a region of interest to examine. In the 2009 data there were 
few high-px objects, so the results here are taken from a mode which searches 
for tracks across the entire tracking detector and is intended for B-physics and 
beam-position determination at L2. 

Offline tracks with jdol < 1-5 mm and \zq\ < 200 mm are matched to L2 
tracks if they are within AR = \/ Aif- + Af/)^ < 0.1. The efficiency is defined 
as the fraction of offline tracks which are matched and is shown in Fig. [16] as a 
function of the track pT- 
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Figure 16: The efHciency for reconstructfoii of a L2 track candidate as a function 
of the pt of the matched offline track in data and Monte Carfo simulation. A 
fit of the threshold curve is superimposed. 



6 Electrons and Photons 

The electron and photon reconstruction and identification algorithms used in 
ATLAS are designed to achieve a large background rejection and a high and 
uniform efficiency over the full acceptance of the detector for transverse ener- 
gies above 20 GeV. Using these algorithms on the 0.9 TeV data, a significant 
number of low-pT electron and photon candidates were reconstructed. The mea- 
surements provide a quantitative test of both the algorithms themselves and the 
reliability of the performance predictions in the transverse energy range from 
the reconstruction threshold of 2.5 GeV to about 10 GeV. 

The electromagnetic calorimeter (EM) consists of the barrel (EMB) and two 
end-caps (EMEC). The barrel covers the pseudorapidity range \ri\ < 1.475; the 
end-cap calorimeters cover 1.375 < \r]\ < 3.2. In the forward direction energy 
measurements for both electromagnetic and hadronic showers are provided by 
the Forward Calorimeter (FCal) in the range 3.1 < \r]\ < 4.9. The hadronic 
calorimetry in the range I77I < 1.7 is provided by the scintillator-tile calorimeter 
(Tile). For 1.5 < \ri\ < 3.2 hadronic showers are measured by the hadronic 
end-caps (HEC), which use LAr with a copper absorber. 

The e/7 algorithms make use of the fine segmentation of the EM calorimeter 
in both the lateral and longitudinal directions of the showers [T] . At high energy, 
most of the EM shower energy is collected in the second layer which has a lateral 
granularity of 0.025 x 0.025 in 77 x space. The first layer consists of finer-grained 
strips in the 77-direction (with a coarser granularity in (jj), which improves 7- 
7r° discrimination. A third layer measures the tails of very highly energetic EM 
showers and helps in rejecting hadron showers. In the range \r]\ < 1.8 these 
three layers are complemented by a presampler layer placed in front with coarse 
granularity to correct for energy lost in the material before the calorimeter. 

The algorithms also make use of the precise track reconstruction provided 
by the inner detector. The TRT also provides substantial discriminating power 
between electrons and pions over a wide energy range (between 1 and 200 GeV). 
The Pixel B-layer provides precision vertexing and significant rejection of photon 
conversions through the requirement of a track with a hit in this layer. 
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6.1 Electron and Photon Reconstruction 

The basic algorithms for electron and photon reconstruction are described in 
detail in Ref. [TB]. The first stage of the search for EM objects is to look for 
significant deposits in the EM calorimeter cells inside a sliding window as it 
is moved across the detector. The size of the sliding window cluster depends 
on the type of candidate (electron, unconverted or converted photon) and the 
location (barrel, end-caps) . The cluster energy is calculated from the amplitudes 
observed in the cells of the three longitudinal layers of the EM calorimeter and 
of the presampler (where present). The calculation sums the weighted energies 
in these compartments, then takes into account several corrections for shower 
depth, lateral and longitudinal leakage, local modulation etc. The weights and 
correction coefficients were parameterized from beam-tests [T] and simulation. 
Electrons are reconstructed from the clusters if there is a suitable match with 
a particle track of pt> 0.5 GeV. The best track is the one with an extrapolation 
closest in (77, 4>) to the cluster barycentre (the energy-weighted mean position) in 
the middle EM calorimeter layer. Similarly, photons are reconstructed from the 
clusters if there is no reconstructed track matched to the cluster (unconverted 
photon candidates) or if there is a reconstructed conversion vertex matched to 
the cluster (converted photon candidates) . "Single track conversions" (identified 
via tracks lacking a hit in the B-layer) are also taken into account. First, 
electron candidates with a cluster \r]\ < 2.47 and photons with cluster |?7| < 2.37 
are selected and investigated (the cluster t] is defined here as the barycentre 
of the cluster cells in the middle layer of the EM calorimeter). Electron and 
photon candidates in the EM calorimeter transition region 1.37 < \r]\ < 1.52 
are not considered. At this stage, 879 electron and 1 694 photon candidates are 
reconstructed in the data with Et above 2.5 GeV. 

6.2 Electron and Photon Identification 

The isolated electron and photon identification algorithms rely on selections 
based on variables which provide good separation between electrons/photons 
and fake signatures from hadronic jets. These variables include information from 
the calorimeter and, in the case of electrons, tracker and combined calorime- 
ter/tracker information. There are three classes of electrons defined: loose, 
medium and tight, and two for photons: loose and tight. The selection criteria 
were optimized in bins of Et and 77, separately for electrons, unconverted and 
converted photons. 

The loose selection criteria are based on the shower shape and are common 
to electrons and photons. For electrons, the medium requirements make use of 
the track information while in the tight ones the particle track selections are 
more stringent and use the particle identification capability of the TRT. For 
photons the tight selection criteria make full use of the EM calorimeter strip 
layer information, mainly to reject merged photon pairs from high energy tt'^'s. 

In the following all reconstructed electron and photon candidates with cluster 
Et > 2.5 GeV at the sliding window level are considered. 
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6.3 Electron Candidates 



Figure [T7] displays, for all of the 879 electron candidates from 384186 events, the 
transverse energy and pseudorapidity spectra. Table |3] presents the percentage 
of these candidates which pass the successive selection criteria both for data 
and simulation. These criteria were not optimized for such low-energy electron 
candidates (see Section 15^ . Both Fig. [T71 and Table [3] show similar behaviour 
in data and simulation. The remaining discrepancies in the first stages of back- 
ground rejection may be related to the small differences observed in shower 
variables (see Section [6.6.11 below) . In Fig. flW b') the drop in efficiency around 
\r]\ = 1.5 corresponds to the barrel/end-cap transition. 



Table 3: The fraction of electron and photon candidates passing the different 
selection criteria, compared to those predicted by Monte Carlo (MC). Statistical 
error are quoted. 



Electron candidates Photon candidates 
Data (%) MC (%) Data (%) MC (%) 



Loose 

Medium 

Tight 



46.5±1.7 
10.6±1.0 
2.3±0.5 



50.9±0.2 
13.1±0.2 
2.4±0.1 



25.4±1.0 

n.a. 
4.1±0.5 



30.5±0.1 

n.a. 
6.6±0.1 



In these figures the Monte Carlo prediction is sub-divided into its two main 
components: hadrons and real electrons. The latter is largely dominated by 
electrons from photon conversions, but also includes a small fraction (^ 3%) 
of electrons from other sources, such as Dalitz decays, and an even smaller 
one (below 1%) of electrons from b, c — >■ e decays. There are twenty electron 
candidates passing the tight selections in the data. Approximately 15% of such 
candidates in the Monte Carlo are from heavy flavour decays. 
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Figure 17: Distribution of cluster Et (a) and |?7| (b) for all selected electron 
candidates. The simulation is normalized to the number of data events. 



6.4 Photon Candidates 

Transverse energy and pseudorapidity spectra for all 1 694 photon candidates 
are displayed in Fig. [THl Table |3] presents the percentage of photon candidates. 
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Figure 18: Cluster Et (a) and \ri\ (b) for all seleeted photon eandidates. The 
simulation is normalized to the number of data events. 



as a function of the selection level applied. Of the selected candidates, 14% are 
reconstructed as converted photons and almost all of these, ~ 98%, are also 
selected as electron candidates. 

The Monte Carlo prediction is sub-divided in this case into four components 
of decreasing importance: approximately 71% of the candidates correspond to 
photons from tt" decay, whereas ~ 14% are from r], i] or oj decays into pho- 
tons; '^ 14% are from other hadrons with complex decay processes and particles 
interacting in the tracker material. At these energies, only a very small fraction, 
~ 1%, of all photon candidates are expected to be primary products of the hard 
scattering. 

6.5 First-level Electron and Photon Trigger Performance 

The LI e/7 selection algorithm searches for narrow, high-iJ^ electromagnetic 
showers and does not separate electrons from photons. The primitives for this 
algorithm are towers which sum the transverse energies of all electromagnetic 
calorimeter cells in A77 x A(/) = 0.1 x 0.1. The trigger examines adjacent pairs of 
towers and tests their total energy against several trigger thresholds. Isolation 
requirements were not yet employed. The lowest threshold e/7 trigger for the 
2009 data-taking period required a transverse energy of at least 4 GeV. 

Clusters consistent with originating from an electron or photon are selected 
by requiring at least 30% of the cluster energy to be deposited in the second layer 
of the electromagnetic calorimeter, where the maximum of an electromagnetic 
shower is expected. The selected e/7 candidates are matched to LI clusters in 
1] and (j) by requiring AR < 0.15. The efhciency is then calculated from the 
fraction of reconstructed clusters which have a matching LI cluster. 

The resulting LI trigger efficiency for the lowest threshold component is 
shown in Fig. 1191 The sharpness of the efficiency turn-on curve around threshold 
agrees with the Monte Carlo expectation. Low energy reconstructed clusters 
occasionally fire the trigger, especially when the coarser granularity used at LI 
merges two separate offline clusters. 
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Figure 19: EfBciency for the lowest threshold LI electromagnetic trigger, a 
nominal 4 GeV, as a function of the uncalibrated offline cluster transverse energy. 
The turn-on is shown for data (solid triangles) and non-diffractive minimum-bias 
simulation (open circles). 
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Figure 20: Fraction of energy deposited by photon candidates with Et > 
2.5 GeV in each layer of the electromagnetic calorimeter for data and simu- 
lation. These fractions are labelled as (a) /o for the presampler layer, (b) /i for 
the strip layer, (c) /2 for the middle layer and (d) /s for the back layer. Frac- 
tions can be negative due to noise fluctuations. The simulation is normalized 
to the number of data events. 
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6.6 Electron and Photon Identification Variables 

6.6.1 Calorimeter Variables 

In this section, various calorimeter-based quantities are displayed for the photon 
candidates. These are preferred to the similar electron distributions because of 
the higher purity. 

Figure [201 illustrates the longitudinal development of the shower in the suc- 
cessive layers of the EM calorimeter, based on the measured layer energies before 
corrections are applied. For the observed photon candidates, which in simulation 
are predominantly from 7r° decays, the energy is deposited in earlier calorimeter 
layers than typical for high energy photons. In the presampler part, the simu- 
lation points are higher than the data for fractions above 0.6. This is at least 
in part because the presampler simulation does not describe the recombination 
of electron-positron pairs by highly ionizing hadrons or nuclear fragments that 
lower the LAr response, an effect which is included in the accordion calorimeter 
simulation. This feature also explains the observed disagreement in the first bins 
for the fractions in the other layers, since the various fractions are correlated. 

Several variables are used to quantify the lateral development of the shower. 
From these, the distribution of W2, the shower width measured in the second 
layer of the EM calorimeter, is shown in Fig. [2lT a). The shower width W2 is 
slightly larger in the data. Preliminary studies show that including the cross- 
talk between neighbouring middle layer cells (^ 0.5%) [T7] in the simulation 
explains part of the observed difference. 

The distribution of two variables used for 77^/7 separation in the tight photon 
selection, i^ratio and Wg^, are shown in Figs. [21T b) and [2ljc). i^ratio is the 
difference of the highest and second highest strip energies, divided by their 
sum. Wss is the shower width measured in three strips around the maximum 
energy strip. For this variable the data show a slightly wider profile than the 
simulation, although in this case the simulation already includes the measured 
cross-talk. In general, all the shower shape variables show good agreement 
between data and simulation. 

6.6.2 Tracking and Track-Cluster Matching Variables 

Electron and converted photon identification rely heavily on tracking perfor- 
mance. Figure [221 illustrates two of the track-calorimeter matching variables 
used in the identification of electron candidates in data and simulation. For 
simulation, hadrons and real electrons are shown separately. Fig. I22r a) shows 
the difference in 77, A771, between the track extrapolated to the strip layer of 
the EM calorimeter and the bary centre of the cell energies in this layer. Fig- 
ure [521[b) shows the difference in azimuth, A(/)2, between the track extrapolated 
to the middle layer of the EM calorimeter and the barycentre of the cell energies 
in this layer. This variable is signed by the charge of the particle to account 
for the position of any radiated photons with respect to the track curvature, 
and an asymmetric cut is applied in the selection. The asymmetric tails at 
large negative values of /S.(j)2 are more pronounced for the electrons than for the 
hadrons. 

Figure [23] shows a comparison of four of the tracking variables between data 
and simulation for all electron candidates. Figures [23|a) and [23l b) show the 
numbers of hits on the tracks in the Pixel and SCT detectors, respectively. The 
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Figure 21: Distributions of calorimeter variables compared between data and 
simulation for all photon candidates with px above 2.5 GeV. Shown are the 
shower width in the middle layer of the EM calorimeter, W2 (a), and the variables 
^-ratio (b) and Wa3 (c), which characterize the shower shape in the first (strips) 
EM layer. The simulation is normalized to the number of data events. 
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Figure 22: Distributions of track-calorimeter matching variables for all electron 
candidates compared between data and simulation, (a) shows the difference in 
T] in the first calorimeter layer (see Section 16.6. 2p and (b) shows the match in 
charge-signed-(^ in the second. The simulation is normalized to the number of 
data events. 
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Figure 23: Distributions of tracking variables for all electron candidates com- 
pared between data and simulation. The number of Pixel (a) and SCT (b) hits 
on the electron tracks are shown, the fraction of high-threshold TRT hits for 
candidates with \ri\ < 2.0 and with a total number of TRT hits larger than 
ten (c), and the transverse impact parameter, do, with respect to the recon- 
structed primary vertex (d). The simulation is normalized to the number of 
data events. 
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fraction of high-threshold TRT hits belonging to the track for electron candi- 
dates with 1 77 1 < 2.0 and with a total number of TRT hits larger than ten is 
shown in Fig. HS^c). At these low energies the transition radiation yield of elec- 
trons is not optimal and yet a very clear difference can be seen between the 
distributions expected for hadrons and for electrons from conversions. Finally, 
Fig. I23f d) shows the distribution of the transverse impact parameter, do, of 
the electron track with respect to the reconstructed primary vertex position in 
the transverse plane; whereas the hadrons in the simulation display a distri- 
bution peaked around zero with a resolution of ~ 100 /im, the electrons from 
conversions have large impact parameters. The agreement between data and 
simulation is good, despite the complications expected at these low energies due 
to material effects and track reconstruction inefficiencies. 
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Figure 24: Ratio, E/p, between cluster energy and particle track momentum (a) 
for electron candidates and (b) for electrons from converted photons. In each 
case candidates withpx above 2.5 GeV in the calorimeter are shown. Sub- figure 
(a) is dominated by real electrons. The simulation is normalized to the number 
of data events. 



Figure IMT a) shows the distribution of the ratio E/p of cluster energy in 
the calorimeter to track momentum for all electron candidates and for data 
and simulation. Electrons from conversions have a broad E/p distribution as 
their shortened tracks have a large momentum error. The hadron component 
peaks at values near unity: this behaviour, due to the selection bias for these 
hadrons, is also observed in the simulation. In a similar fashion. Fig. IMT b') 
shows the E/p ratio of the reconstructed converted photon candidates, where 
the converted photon momentum is estimated from the combination of the par- 
ticle momenta for double-track conversions and from the particle momentum 
measurement available for single-track conversions. Approximately 20% of the 
converted photon candidates are reconstructed as single-track conversions in 
this kinematic regime. Both the electron dominated and the hadron dominated 
distributions show good agrement with the simulation. 

6.6.3 Use of the TRT for Electron Identification 

As already discussed in Section [^31 the electron candidate data sample is ex- 
pected to consist predominantly of two components: charged hadrons misre- 
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constructed as electrons and electrons from photon conversions. These two 
components can be separated by using the measured fraction of high-threshold 
TRT hits on the electron tracks (see Fig. B^ c)). To perform such a measure- 
ment, the electron candidates are required to lie within the TRT acceptance, 
i.e. \ri\ < 2.0, and to have a reconstructed track with a total of at least ten TRT 
hits. 

The distribution of the fraction of high threshold hits has been fitted in 
20 bins between and 0.5 to extract the number of hadrons and electrons 
observed in the data. This relies upon the modelling of the response of the TRT 
to electrons and pions in the simulation. The sample of electron candidates 
considered here is predicted to contain 494±26 electron candidates which are 
actually hadronic fakes and 226±21 genuine electrons. 

Two examples of comparisons between the shapes of variables extracted 
for each of the two components, using the method described above (on each bin 
individually), and the shapes predicted for each component are shown in Figs. 1251 
and [26l respectively, for two of the most sensitive variables: the fraction of 
the cluster energy measured in the strip layer and the ratio E/p. The E/p 
distribution for electrons in both data and simulation shows a peak close to unity 
and a tail at large values from bremsstrahlung losses in the tracker material. The 
error estimates in these de-convolved plots come from toy Monte Carlo trials 
and their size reflects the power of the TRT detector for electron identification. 
The rates of electrons and hadrons and the relevant distributions agree with the 
Monte Carlo simulation for each species illustrating the quality of the simulation 
modelling. 
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Figure 25; Distribution of the energy fraction in the strip layer of the 
EM calorimeter as extracted from data compared to the truth from simula- 
tion. The results are shown for both components of the electron candidates: 
electrons from conversions (a) and hadrons (b). The simulation is normalized 
to the number of data events. 

6.7 Photon Conversions 

An accurate and high-granularity map of the inner detector material is nec- 
essary for a precise reconstruction of high-energy photons and electrons. The 
location of the conversion vertex can be used as a tool to map the position and 
amount of material of the inner detector. In the following, photon conversions 
are selected using only information from the inner detector, enabling the use of 
very low momentum particle track pairs. In addition, conversions give a source 
of electrons from which the TRT detector response can be determined. 
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Figure 26: Distribution of the E/p as extracted from data compared to the truth 
from simulation. The results are shown for both components of the electron 
candidates: electrons from conversions (a) and hadrons (b). The simulation is 
normalized to the number of data events. 



The conversion reconstruction algorithm is described in detail elsewhere [16j . 
In the following, the basic steps of the algorithm are recalled together with 
an updated list of selection criteria. The algorithm begins by selecting single 
particle tracks with transverse momentum pt> 500 MeV. These tracks must 
have a probability of being an electron of more than 10%, calculated using the 
particle identification capability of the TRT, see Section [531 

Conversion candidates are then created by pairing oppositely charged par- 
ticle tracks. The tracks are further required to be close in space and to have 
a small opening angle. The selected particle track pairs are then fitted to a 
common vertex with the constraint that they be parallel at the vertex. The 
final set of conversion candidates is selected based on the quality of the vertex 
fit which must have x^ smaller than 50. 

The tracks used for the reconstruction of conversions may be stand-alone 
TRT tracks, or they may include silicon hits. In the data, 3 662 vertices, 6.7% 
of the total, have two tracks with silicon information, to be compared with 
10.4% in the simulation. This class of vertices have much less background than 
the total and the following results are drawn from them. Some properties of 
the candidates in data and Monte Carlo simulation are shown in Fig. [571 Given 
the complexity of the reconstruction of converted photons and the impact of 
bremsstrahlung of the electrons in the tracker material, the consistency between 
the data and the simulation for the selection variables is good. 

To measure the inner detector material the selection requirements are tight- 
ened to >90% TR electron probability and vertex x^ < 5. 

Figure [^5] shows the location in radius and rj of conversion vertices. The 
simulation was normalized to the same number of conversions as in the data 
and the agreement in shape is satisfactory. 

The amount of material, in multiples of the radiation length Xg, that the 
photons traverse can be calculated from the fraction of photon conversions seen 
in it given the reconstruction efficiency. The combinatorial background and the 
error in determining the conversion radius must be accounted for. To remove 
the dependence on the absolute flux of photons and the overall reconstruction 
efficiency, the rate is normalized to that seen in a well-understood reference ma- 
terial volume, which is chosen to be the beam pipe. As shown in Fig. [28l there 
were only 9 conversions in this reference volume and so the absolute material 
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Figure 27: Comparison between converted photon candidates, for which both 
tracks have sihcon hits, in data and non-diffractive minimum-bias Monte 
Carlo simulation, (a) Opening angle in the rz plane between the two tracks 
(A(l/tan6')); (b) 3D distance of closest approach between the two tracks, dr. 
The distributions are normalized to the same number of conversion candidates 
in data and Monte Carlo simulation. 

determination has errors of at least 30%. The agreement between data and 
Monte Carlo is presented in Table |H 



6.8 Reconstruction of 7r° and rj Mesons 

For the analysis presented in this section, cells from the four layers are com- 
bined to form a cluster of size A77 x A0 = 0.075 x 0.125, which corresponds 
to an area of 3 x 5 cells in the middle layer of the EM calorimeter. The EM 
cell clusters are reconstructed with a seed cell threshold |£'ccii| — 4cr (where a 
corresponds to the electronic noise in the cell) and with a cluster transverse 
energy Et > 300 MeV [K]- These clusters are used as photon candidates for 
tt" and 77 reconstruction. 

The standard parameterization of energy response discussed in Section 16.11 
was performed for photons with Et > 5 GeV. For the present study a dedicated 
parameterization was extracted from the minimum-bias simulation sample using 
low-energy photons coming only from tt'^s. 



6.8.1 Extraction of tt° -^ 77 Signal 

In order to extract the tt'^ signal from the combinatorial background, well mea- 
sured photons were selected inside an acceptance of |r/| < 2.37, excluding a 
transition region 1.37 < |ry| < 1.52. The fraction of energy in the first layer. 
El /{El + E2 + E3), was required to be larger than 0.1 and the clusters were 
required to have a transverse energy, Et , above 400 MeV. 

All pairs of photons with p^^^ > 900 MeV are selected. There are about 
8 X 10^ of these in the data. 
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Figure 28: Distribution of conversion candidate radius, (a), and 77, (b). The 
points show the distribution for data; the open histograms, the total from the 
Monte Carlo simulation and the filled component shows the expected contri- 
bution of true photon conversions. The contribution from the Dalitz decays of 
neutral mesons is shown in sub-figure (a). The Monte Carlo simulation is nor- 
malized to number of conversion candidates in the data, although in subsequent 
analysis normalization is to the number in the beam pipe. 

6.8.2 7r° Mass Fit 

The invariant mass distribution of the photon pairs is shown in Fig. [29] for 
both data and Monte Carlo. The diphoton mass distribution is fitted using a 
maximum-likelihood fit. The signal is described by the sum of a Gaussian and 
a "Crystal-Ball function" [TS], which are required to have the same mean. The 
combinatorial background is described with a 4"^ order Chebyshev polynomial. 
The parameters of the signal and the background normalization are varied in 
the fit to the data, while the parameters of the polynomial were extracted from 
the Monte Carlo. 
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Figure 29: (a) Diphoton invariant mass distribution for the 7r° selection for 
data and Monte Carlo. The Monte Carlo is normalized to the same number of 
entries as the data, (b) Invariant mass distribution from one converted and one 
imconverted photon. The data are represented by points and the Monte Carlo 
simulations are shown as histograms. 
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Table 4: A^roco is the number of reconstructed conversions in each layer, and 
^/^odata ^-iid X/Xqmc represent the amount of material in the different 
volumes estimated from data and Monte Carlo, normalized by the number 
of reconstructed converted photons in the beam pipe, whose material is 
assumed to be correct. The normalization introduces an additional statistical 
uncertainty of 30% on X/Xo^ata- 
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The fitted n° mass is 134.0±0.8 MeV for the data and 132.9±0.2 MeV for 
the Monte Carlo where the errors are statistical only. The mass resolution in 
the data is 24.0 MeV, to be compared with 25.2 MeV in the simulation, and the 
number of tt'^'s is (1.34 ± 0.02) x 10^. This fit is sensitive to the modelling of 
the background shape near the tt'^ mass. Varying the background shape under 
the peak leads to a differences of up to 1% for the fitted 7r° mass, up to 10% 
for the fitted tt'^ mass resolution and up to 20% for the fitted total number of 
signal events. 

The 1% agreement of energy scale between data and Monte Carlo is well 
within the 2 - 3% uncertainty on the energy scale transported from test-beam 
data analysis. The 1.5% discrepancy of the mass found in Monte Carlo with 
respect to the PDG nominal 7r° mass is consistent with the accuracy (as evalu- 
ated with simulation) of the cluster calibration procedure used for the low-energy 
photons, and the 1% uncertainty arising from the background modelling. 

The converted photons reconstructed in Section [^?71 can also be used to search 
for the TT*^. This is done here using one photon reconstructed in the calorime- 
ter, with a track veto applied and one conversion candidate. The conversion 
candidates are required to have four silicon hits on both tracks and must be in 
the same hemisphere as the calorimeter cluster. Figure l^^f b) shows the je~^e~ 
mass spectrum; the tt'^ peak is clearly visible. 

The uniformity of the EM calorimeter response was studied in ten i] bins, 
where both photons are in the same bin. The diphoton mass distribution in 
each ry bin is fitted separately with the background shape constrained from 
simulation. The reconstructed 7r° mass is constant within 3% for both data and 
Monte Carlo for all rj bins, and the ratio of data to Monte Carlo is consistent 
within the 2% statistical uncertainties. 

6.8.3 Extraction of the rj — >■ 77 Signal 

The number of ry ^- 77 events is expected to be one order of magnitude smaller 
than tt'^ —>■ -f-f in the minimum-bias event sample. Therefore, the combinatorial 
background contribution in the 77 mass region needs to be significantly reduced. 
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Figure 30: Diphoton invariant mass spectrum with tighter selection criteria to 
extract the 77 peak with the fit superimposed to the data. The Monte Carlo 
simulation sample is normalized to the number of entries in the distribution for 
data. 



This can be achieved by adding the following criteria to the tt" analysis: 

• Tighter kinematic selections: E^'""*'^'' > 800 MeV, pPf'' > 2200 MeV. 

• A track veto: no track, extrapolated into the calorimeter, should be within 
-0.1 < ((/)cius - 4>cxtr) < 0.05 and jrycius - Jycxtrl < 0.05 of the cluster being 
considered. 

The diphoton invariant mass spectrum of this sample is shown in Fig. [301 for 
both data and Monte Carlo. In addition to the 7r° peak, the 77 — >■ 77 signal can 
be observed on top of the combinatorial background. The mass spectrum was 
fitted using the sum of a Gaussian and a Crystal-Ball function with the same 
mean for the 7r° peak, a Gaussian for the rj peak and a 4*^ order Chebyshev 
polynomial for the background. The tt*' and background shape parameters are 
taken from the Monte Carlo simulation while their normalizations are free in 
the fit, as are the parameters of the Gaussian describing the r/ peak. 

As can be seen from Fig. 1301 the number of 77 candidates per photon pair 
agrees between the data and the Monte Carlo simulation. The 77 mass extracted 
from the data, 527 ±11 (stat) MeV, agrees with the mass obtained using the 
same fitting function on the Monte Carlo simulation, 544 ± 3 (stat) MeV, within 
the statistical and energy scale uncertainties. 



7 Jets 

Many of the final states which will be studied in high energy collisions contain 
jets of hadrons produced by strong interactions. The ATLAS analysis chain 
applies the same jet algorithm to the 0.9 TeV and 2.36 TeV collision data and to 
the Monte Carlo simulation. The following comparison between data and Monte 
Carlo simulations should not be taken as a precise analysis of the underlying 
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physics in the simulation, but rather as an assessment of the general behaviour 
of the detector and software chain (reconstruction and simulation) . 

Results are presented using clusters of calorimeter cells calibrated to cor- 
rectly measure the energy deposited by electrons and photons in the calorime- 
ter. This is known as the electromagnetic scale. There was, at this stage, no 
allowance for energy loss in inert material. From this starting point jets were 
reconstructed using the anti-ky algorithm [12], which is safe against infrared 
and collinear divergences. The parameter R, which controls the size of jets in 
the rj — (p plane, was set to R=0.6. 

7.1 Jets from Calorimeter Clusters 

The inputs to the jet algorithm are topological clusters [16] which attempt to 
reconstruct the three-dimensional shower topology of each particle. These clus- 
ters were built starting from seed cells with energies |i?i| > ianoise, where Cnoise 
is the electronic noise measured by iteratively gathering neighbouring cells with 
\Ej\ > 2anoise and, in a final step, adding all direct neighbours of these accumu- 
lated secondary cells. The noise in the EM calorimeter, for example, was in the 
range of 10-40 MeV per cell, depending upon the compartment and pseudorapid- 
ity. Approximately 0.1% of all cells were classified as noisy and were removed. 
Clusters built from the remaining cells were then used to create jets, which have 
to satisfy: p!^ >7 GeV and \ri\ <2.6 where p^r^ is the transverse jet momentum at 
the electromagnetic scale. In order to remove cosmic muons and some residual 
effects from cells in the calorimeter that exhibit large noise fluctuations the jets 
are required to pass quality criteria. Furthermore, if the jet energy corrections 
compensating for excluded calorimeter regions exceed 20%, the corresponding 
candidates are not considered. Figure |3I] presents example distributions of the 
internal structure of the jets, namely the number of topological clusters, and 
the fraction of the jet energy carried by each of them. 
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Figure 31: Distributions of (a) the number of clusters per jet and (b) the fraction 
of energy per cluster for jets reconstructed with topological clusters using the 
anti-ky algorithm with R=0.6. 

Figure [35] shows the jet p^ in data and Monte Carlo simulation, normalized 
to the number of jets in data. Figure|32jb) shows the difference of the azimuthal 
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Figure 32: Distributions of (a) px for all jets and (b) Acf) for events with two 
or more jets. These are shown for jets reconstructed with topological clusters 
using the anti-k^ algorithm with R=0.6. 

angle of the two leading jets (A0) in events with at least two reconstructed jets 
in data. The distribution peaks at tt, corresponding to a topology where the two 
reconstructed jets are back-to-back in the transverse plane. In all distributions, 
the agreement between data and simulations is good, demonstrating that the 
description of the material and detector response in the simulation provides an 
adequate model of basic jet quantities. 

7.2 Performance of the First Level Jet Trigger 

The efficiency of the LI jet trigger has been studied. Jets were reconstructed 
as before, but with a lower energy threshold, E^^ > 4 GeV. For the efficiency 
determination the jet signatures identified by the first level calorimeter trig- 
ger were matched to those of reconstructed jets by requiring AR < 0.6. LI 
jet signatures are based on sums of trigger towers in the electromagnetic and 
hadronic calorimeter both calibrated to the electromagnetic scale. Trigger tow- 
ers are formed by analogue summation on the detector and mostly have a size of 
0.1 X 0.1 in 77 X (/), being larger in parts of the end-caps and in the FCal. For the 
lowest threshold jet trigger the trigger towers are summed over 0.4 x 0.4 in 77 x (/> 
and all higher threshold jet triggers use 0.8 x 0.8 in r] x cj). The special treat- 
ment of the lowest threshold jet trigger is motivated by its higher susceptibility 
to noise. 

The trigger efficiency is calculated as the fraction of reconstructed jets pass- 
ing the quality requirements described which have a matched trigger jet. The 
results for the two jet triggers with 50% efficiency at around 15 GeV and 20 
GeV are shown for data and Monte Carlo simulation as a function of the re- 
constructed jet transverse energy, E';^* in Fig. [321 where i?;^* is given at the 
electromagnetic scale. The curves are fits to the simulation using a standard 
trigger turn-on parameterization with an error function. The data are in agree- 
ment with the fit, which shows that the thresholds are as they are expected to 
be. This is a step to understanding the initial jet selection performance of the 
first-level calorimeter trigger. 

7.3 Tau Studies Using Jets 

Jets from QCD processes form the largest background for the reconstruction 
of hadronically decaying tau candidates. Therefore, even though the actual 
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Figure 33: LI jet trigger efficiency for the triggers with 50% efficiency at around 
15 GeV and 20 GeV for data (soHd) and simulation (open) together with a fit 
to the Monte Carfo using a standard threshold function (see text). 
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Figure 34: (a) The electromagnetic radius Rem (see text) of the inclusive 
reconstructed tau candidates, (b) The same variable with a tightened selection 
requiring di-jet events. The Monte Carlo is normalized to the same number of 
candidates as in the data. 
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number of tau leptons in the 2009 data is expected to be negligible, basic track 
and cluster distributions have been studied using jets with emphasis on the 
variables of importance for tau reconstruction and identification. In Fig. 1341 
the electromagnetic radius, Rem, (the energy- weighted mean AR of the jet 
components from the seed cell) is shown for all tau candidates with E^r^ >7GeV 
on the left and for a subset with a tight di-jet selection, which are more likely 
to be taken as r candidates, on the right. Good agreement between data and 
simulation is observed, which gives confidence in the results obtained from earlier 
performance estimates based on Monte Carlo simulations. 

8 Missing Transverse Energy 

A reliable measurement of the missing transverse energy, EJ^^^'^ , is a key ingredi- 
ent for many important analyses. This study considers E^^^^ reconstructed from 
calorimeter information only. As this measurement involves summing calorime- 
ter cells over the whole detector, it is sensitive to detector and reconstruction 
effects. In particular events with rare unexpected noise contributions tend to 
appear in the tail of the i?™'^^ spectrum. 

^miss reconstructed with the calorimeter is derived from the vector sum 
of the transverse energies of the cells. Because of the high granularity of the 
calorimeter (about 187000 cells), it is crucial to suppress noise contributions 
to E^^^^, i.e. to limit the number of cells used in the sum. This is done by 
only using cells belonging to three-dimensional topological clusters defined in 
Section 17.11 About 800 and 2 500 cells on average are included in such clusters 
in randomly triggered and collision events, respectively. The 0.1% of the cells 
classified as noisy are removed. 

The sensitivity to noise can be best studied in randomly triggered events, 
where minimal energy is deposited in the calorimeters. The E^^^^ distribution 
of these events is shown in Fig. 1351 demonstrating the level of tails in randomly 
triggered events. 
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Figure 35: Distribution of E^^^^ as measured in data from randomly triggered 
events. Only cells belonging to topological clusters are included in the calcula- 
tion; their energies are calibrated at the EM scale. 
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In soft proton-proton collisions, no true Ei^^^^ is expected. This is confirmed 
by the Monte Carlo simulation. Unlike in randomly triggered events, total 
transverse energies (^ Et) up to 100 GeV are deposited in the calorimeter for 
minimum-bias events in the present data set. Figure [35] shows the measured 
E'^''''' distributions as an example. For both E'^''''' and £'^'"", the RMS values 
are about 1.4 GeV and 1.8 GeV for 0.9 TeV and 2.36 TeV centre-of-mass energies, 
respectively. These values are higher than in randomly triggered events because 
the finite resolution in the presence of real energy contributes to the width. 
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Figure 36: Distribution of E"^^^^ (a,b) and E^^^^ (c,d) as measured in data from 
minimum-bias events (dots) at 0.9 TeV (a,c) and 2.36 TeV (b,d) centre-of-mass 
energies. In the calculation only topological cluster cells are used, with energies 
calibrated at the EM scale. The expectations from Monte Carlo simulation are 
superimposed (histograms) and normalized to the number of events in data. 

The E^^^^ distribution, also shown in Fig. [36l is found to be satisfactory at 
this early stage. There are no events outside the borders of the figure. The two 
data events with large EJ^^^'^ are consistent with arising from out-of-time energy 
in the detector. At least one of these appears to be a cosmic ray event. Such 
events are not included in the Monte Carlo simulation sample. 

A more quantitative evaluation of the E^^^^ performance can be obtained 
from a study of the E^^^^ and E^^^^ resolutions as a function of the total 
transverse energy ^ E^ in the event. The resolutions are expected to be pro- 
portional to ^/^ET■ The resolutions observed in the ATLAS data at both 
centre-of-mass energies are presented as a function of \/^Et in Fig. [ST] A 
very good agreement between data and Monte Carlo simulation is obtained at 
both centre-of-mass energies. The E!^^^^ resolution can be parameterized as 
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^(^miss^£;miss-)^ 0.37 X y/J^^T at the EM scale!, with a negUgible statistical 



error. 
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Figure 37: E!f^^^ resolution as a function of the total transverse energy (^ Et) 
for minimum-bias events. The line shows a fit to the resolution expected from 
the Monte Carlo simulation and the full dots (open squares) represents the 
resuhs with data at 0.9 (2.36) TeV. £;^*"", E^'^", J2 Et are computed with 
topological cluster cells at EM scale. 



9 Muons 

The calorimeters are surrounded by the muon spectrometer j5], which was de- 
signed to provide a trigger and accurate standalone reconstruction for muons 
with pt from several GcV up to a few TeV. In contrast, the muons studied in 
the sample collected at y/s—0.9 TeV are of relatively low px- The air-core toroid 
system, with a long barrel and two inserted end-cap magnets, generates an av- 
erage field of 0.5 T (1 T), in the barrel (end-caps) over a large volume. Multiple 
scattering effects are thereby minimized, and excellent muon momentum reso- 
lution is achieved with three layers of high precision tracking chambers. Over 
most of the 77 range, tracks are measured by Monitored Drift Tubes (MDT). For 
2 < I77I < 2.7, Cathode Strip Chambers (CSC) with higher granularity are used 
in the innermost layer, to withstand the demanding rate and background condi- 
tions. In addition, the muon spectrometer includes dedicated trigger chambers 
with nominal timing resolutions between 1.5 and 4 ns. They are composed of 
Resistive Plate Chambers (RFC) in the barrel and Thin-Gap Chambers (TGC) 
in the end-cap regions. Besides providing trigger signals they also measure 
the muon coordinate in the direction orthogonal to that determined by the 
precision-tracking chambers. 

The muon analysis uses a somewhat smaller data set than other analyses 
as the toroid system was not always operational. In addition, it was checked 
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Figure 38: The distribution of 77 and p of reconstructed muons in the 0.9 TeV 
data. The simulation distributions are normahzed to the number of entries in 
data. 

that the MDTs, TGCs and RPCs were all operating normally. The algorithms 
used for muon reconstruction combine tracks from the muon systems and the 
inner detector, and are developed from those described in Ref. [TC]. For the 
results presented here, two independent algorithms are used and only candidates 
selected by both are accepted. This selects a total of 50 muon candidates. Raw 
kinematic distributions for these candidates are presented in Fig. 1381 The muon 
spectrum observed is soft and strongly peaked in the forward direction, where 
the momentum of the muons more often exceeds the 3.2 GeV needed to penetrate 
through the forward calorimeter. The kinematic distributions are compared to 
the predictions from minimum-bias simulation with the Monte Carlo normalized 
to the number of muons found in data. Within the large statistical uncertainties 
good agreement is found indicating a reasonable understanding of the initial 
performance of the ATLAS muon spectrometer. 

The muon trigger, designed to select high-px muons, has a limited acceptance 
for the muon tracks reconstructed offline in 2009. Of the 38 muons in the end- 
cap regions, 13 were triggered at LI by the TGC. Of the 12 muons in the barrel, 
one was triggered at LI by the RFC with the correct timing, another 9 were 
triggered with -j-l bunch crossing misalignment during the timing adjustment 
phase, while the other 2 muons were outside the trigger acceptance. Only one 
muon passed the full trigger chain up to the EF combined trigger after applying 
the pt > 4 GeV selection. The muon momenta and directions measured by the 
L2 and EF are in good agreement with the offline measurement. 

10 Conclusions 

The overall performance of the ATLAS detector at the LHC was established in 
first collision data at centre-of-mass energies of 0.9 TeV and 2.36 TeV. Although 
the detector has not been optimized for the low energy particles studied, its 
performance was found to be remarkably good, particularly in view of the early 
stage of data taking. 
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The overall data-taking efficiency was about 90% and the subdetectors were 
typically 99% operational. The entire computing infrastructure of Trigger/DAQ 
was immediately effective. Collision candidates were selected with a negligible 
background level. 

The tracking detectors and electromagnetic calorimeters have, by the na- 
ture of the data set, been the most extensively tested components and they 
perform well. The hit efficiencies, resolutions and particle identification capa- 
bilities of the tracking detector are well modelled by Monte Carlo simulations. 
The discrimination between electrons and pions using transition radiation was 
demonstrated and outperforms any previous colliding beam experiment. The 
results of various studies suggest that the material distribution in the inner part 
of the tracker is well understood, at a level of a few percent of the total in the 
barrel silicon systems. 

The momentum scale linked to the Kg mass is known at the per mille level 
and the calibration of the electromagnetic calorimetry in the region of the tt^ 
mass was checked at the level of 1%. Electron and photon reconstruction was 
extensively tested and performs well. The candidates typically have transverse 
momenta well below those for which reconstruction and identification algorithms 
were optimized, but their properties are shown to be in good agreement with 
simulation. Good calorimeter performance was also demonstrated by the mea- 
surement of the resolution of the missing transverse energy, which follows closely 
the expectations from Monte Carlo over the entire energy range probed. 

The muon system was not extensively tested with this data set, but performs 
as expected within the precision available. It was well tested in cosmic ray data 
taken in 2008 g]. 
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R. Konoplich^"^, S.P. Konovalov^"*, N. Konstantinidis^"^, A. Kootz"^^ 

S. Koperny^'^, S.V. Kopikov^^s^ ^ Korcyl^*, K. Kordas^^^, V. Koreshev^^®, 
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C. Kourkoumelis^, A. Koutsman^"'"', R. Kowalewski-'^^^, H. Kowalski'*^, 

T.Z. Kowalski37, W. Kozaneckii^e^ A.S. Kozhini^s, V. Kral^^^^ 

V.A. Kraniarenko^^, G. Kramberger^^, O. Krasel^^, M.W. Krasny^*, 
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L. Nicolas^^Q^ q Nicoletti"'^, B. Nicqucvcrt^s, F. Niedercorn"^ M. Magic's.*, 
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H. Okawai62^ Y. Okumura^°\ T. Okuyamai^", M. Olcese^O'^, 

A.G. Olchevski^s, M. Oliveirai^*^^ d. Oliveira Damazio^^, C. Oliver^o, 

J. Oliver^"^, E. Oliver Garcia^^e, D. Olivitoi^o, M. Olivo Gomez^^, 

A. Olszewski's, J. Olszowska's, C. Omachi^^^", A. Onea^^, A. Onoire^^*'', 
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O.K. Oye^^^ v.E. Ozcan'^'^, K. Ozone^'^, N. Ozturk^, A. Pacheco Pages", 

C. Padilla Aranda", E. Paganisms, C. Pahl^^, F. Paige^^, K. Pajchel"^, 

S. Palestini29, J. Palla^^, D. Pallin^^^ ^ Palniai24a^ j_d. Palmeri^, 
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C. Puigdengoles", J. Purdham*^ M. Purohit^*'^, P. Puzo"^ 
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